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A note on the power divergence in lattice calculations of 

AI = 1/2 K -> tttt amplitudes at M K = M n 
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Abstract: 



C*") ■ In this note, we clarify a point concerning the power divergence in lattice calculations of AI = 1/2 
K — > 7T7T decay amplitudes. There have been worries that this divergence might show up in the Minkowski 
amplitudes at Mk = with all the mesons at rest. Here we demonstrate, via an explicit calculation 
in leading-order Chiral Perturbation Theory, that the power divergence is absent at the above kinematic 
point, as predicted by CPS symmetry. 



PACS numbers: 11.15.Ha,12.38.Gc,12.15Ff 



The subtraction of a power divergence, which arises via the mixing of dimension-six four-fcrmion 
operators with those of lower dimension has been one of the central issues in lattice calculations of 
AI = 1/2 K — > tttt amplitudes. This power divergence is of course unphysical, and can be related to a 
shift of the vacuum due to the inclusion of the weak interaction in Chiral Perturbation Theory (%PT) 
Q] [21 [31 S| - It results in the so-called tadpole operators, which contribute to the processes K° — > |0) and 
K° — > |0), in xPT with weak interactions. 

As argued in Ref. [S], this power divergence should be absent for K — > tttt amplitudes when m s = 
md = rn u {m u> d,s are the masses of u, d and s quarks), due to the exact CPS symmetry yQ of the four- 
fcrmion operators that mediate K — > tttt decays. In Ref. [H], it was argued that the power divergence 
indeed does disappear in Euclidean space at Mk = M^. However, a naive calculation in Minkowski space 
suggests that this power divergence might still be present at Mk = M v when all mesons are at rest. The 
issue is relevant, as it has been proposed that this unphysical kinematic point can be used to extract the 
low-energy constants relevant for AI = 1/2 K — ► tttt to order p 4 in xPT [7HH]. 1 

In this note, we show, via an explicit calculation in xPT, that also in Minkowski space the power 
divergence is not present in AI = 1/2 K — > tttt amplitudes at Mk = M w , with all mesons at rest. Since 
it has already been argued in Ref. ||| that the AI = 1/2 K — > tttt amplitudes in partially quenched xPT 
at the kinematic point Mk = M„ suffer from problems related to the lack of unitarity [51 1101 ITT] , we 
concentrate here on full QCD. Our conclusions on the power divergence will however not change in the 
(partially) quenched case. 

To simplify the discussion, we only consider weak operators in the (8, 1) irrep of SU(3)l X SU(S)r. 
The weak mass operator in this irrep at 0(p 2 ) in the chiral expansion is 

0f 1] = a 2 x {2B Tr [A 6 (A^S + ^ M)} } , (1) 

where a 2 is the (power-divergent) low-energy constant associated with this operator, B = — (0\uu + 
dd\0) / f 2 (in the chiral limit), X & is a Gell-Mann matrix, Ai is the quark-mass matrix and S is the 
standard non-linear Goldstone field. 

We first observe that CPS symmetry implies that the parity-odd part of this operator is proportional 
to m s — m-d- In fact, 

Of' 1] = a 2 x {B a (m s + m d )Ti [A 6 (E + S f )] + iB Q (m s - m d )Tr [A 7 (E - E f )] } . (2) 

Therefore, at m s = rrid the parity-odd part of the operator vanishes, and thus its K — > 7T7t matrix element 
should vanish as well for Mk = M w . This was confirmed by an explicit calculation in Euclidean space 
(as reported in Ref. [H]), and should be true in Minkowski space as well. 



At leading order in the chiral expansion, 2 ' contributes to the K — * tttt amplitudes via the 
diagrams in Fig. ^ where the gray circles represent the weak mass operator, and the square is from the 
leading-order strong chiral Lagrangian. In diagram (b), there is a pole associated with the K° propagator 
(the dashed line in Fig. [IJ, which takes the form 



when all the other three on-shell particles are at rest. For fixed Mk ^ M V: one may take e — > at any 
stage of the calculation, since the denominator of Eq. J2J does not vanish in that case. However, for 
Mr = Mm the it prescription is needed in order to define the propagator, and should be taken to zero 
only at the end of the calculation. In that case, one finds: 

(tt+tt-IO^IK^m^ = lim j M Hm^ a 2 x (M| - M 2 ) 

1 It follows from our analysis that the low-energy constants c*2 and 2 5 should n °t appear in Eq. (31) of Ref. 
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Figure 1: Diagrams involving the weak mass operator at the lowest order in the chiral expansion for the 
AI =1/2 K — > 7T7T amplitudes. The gray circles represent the operator 0\ ' , and the square is the 
K°K° — > 7r + 7r~ vertex from the lowest-order strong chiral Lagrangian. The dashed line in (b) indicates 
that the K° could be off-shell, while all the other mesons are always on-shell. 

which indicates that there is no need to perform the subtraction of a power divergence. 

Let us discuss this claim in some more detail. We begin by noting that it was shown long ago ^ dj 
that if the weak mass term J d^xO^' 1 ^ is treated as a perturbation to the strong chiral Lagrangian, it 
does not have any observable effect. However, here we consider the unphysical situation of an energy 
non-conserving matrix clement of 2 (corresponding to the insertion of this operator at a fixed time, 
see below), and the above consideration does not apply. 

The factor (M K — M%) on the right-hand side of Eq. Q originates from the CPS symmetry of the 
operator (c.f. Eq. ©), while the quantity in the square brackets is determined by the kinematics. This 
latter quantity indeed diverges in the limit — > M„ (and e — > 0). That this is exactly what one expects 
to happen follows because the K° propagator in Fig. ^b) goes on-shell without being amputated. In 
fact, Fig. ^b) also represents the process of K°-K° scattering into ir + -ir~, but in that case in order to 
obtain a finite amplitude, the LSZ reduction formula tells us to amputate the K° external leg, before 
putting it on-shell. Since in our case this leg is not amputated, the diagram is divergent in the on-shell 
limit. In the case that K°, ir + and ir~ are all at rest, this u K a on-shell" point coincides with the limit 
Mk — ► and CPS symmetry prevents the divergence from happening: the amplitude actually vanishes 
at M K = M w . 

However, one may consider the following more general situation. Consider for instance kinematics 
with K° at rest but ir + and ir~ carrying spatial momenta p and —p respectively. In that case, the "K° 
on-shell" point is at Mk = -EV = + \p\ 2 , and diagram ^b) is proportional to p 1 jt at this point. 

The extra factor (M K — M%) clearly does not help in this case, and the divergence occurs of course for 
exactly the same reason as described above. 

We gain more insight by considering the amplitude in position space, as in Fig. Gib). This diagram 
contains a factor Q~ lM K\t^-t s \ f rom the K° propagator, where t w is the location (in time) of the weak 
operator 0\ ' (taken as t w = in the diagram), and t s is the location of the strong vertex. The LSZ 
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Figure 2: Diagrams involving the weak mass operator at the lowest order in the chiral expansion for 
the correlator (OlTr+Tr - ^^ 8 ' 1 ^ ^). The gray circles represent the weak mass operator 0^' X \ and the 
square is the K°K° — > tt + it~ vertex from the lowest-order strong chiral Lagrangian. The dashed line in 
(b) means K° could be off-shell, while all the other mesons are always on-shell. The weak operator is at 
the space-time origin. K is created at (jk\ and the pions are annihilated at t v (t w ) in Minkowski 
(Euclidean) space. 



prescription for this K° line corresponds to taking a Fourier transform with respect to t w , and putting 
the corresponding momentum on-shcll. For this to work, the integral over t w needs to be regulated 
by replacing Mk — ► Mk — «e, and this is precisely what leads to the it prescription in Eq. Q. It 
follows that the divergence encountered here is regulated by considering the amplitude at finite t w (by 
time-translation invariance we may choose t w = 0). This is of course what one does anyway in a lattice 
computation of this amplitude. It is therefore instructive to consider this amplitude in position space 
rather than momentum space [5], which is what we will do next. 



Since we take all our mesons to have vanishing spatial momentum, we will consider the relevant 
correlators in the time-momentum representation, i.e. study the correlators as functions of 3-momentum 
and time. In this setup, a free meson propagator with energy Ep = ^Jm? + \p\ 2 (m is the mass and p is 
the 3-momentum of the meson) is 



-iEp\t\ 

-%|t| 



2Ep 



(Minkowski) 



(Euclidean) 



where t (t) is the Minkowski (Euclidean) time. The time dependence of the Minkowski expression is of 
course in accordance with the ie prescription of Eq. 0. We now consider the Minkowski correlator, 

C 2 = (0| 7 r+(t 7r )7r g -( iw )of 1) (0) J Pr 9(^)|0), (5) 
and its Euclidean counterpart, 

C 2 = (OWMt^MO^^O^tkM. (6) 
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For simplicity, wc choose to annihilate the two pions at the same time, and assume that Ik{tk) < 
and ^(tv) > 0. The weak operator is inserted at time t w = (t w = 0). All particles are at rest (as 
indicated by the subscripts 0). The relevant diagrams for the above correlators are shown in Fig. [21 In 
the following, we only present the result in Minkowski space, but stress that the calculation in Euclidean 
space is virtually identical, and leads to the same conclusion UJ. 

The contribution from diagram Ofa) to the correlator C2 is: 
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(2M K )(2M„)(2M„) 



(7) 



while diagram |2Ib) leads to 
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{Ml - Ml) 



Aict2 

W V "~ K (2M K )(2M K )(2M 7r )(2M n ) 

dt e - iM K\is-t K \ e -iM K \t 3 \ c -2iM„\t w -t 3 1 



x {M 2 K [1 + e(t s )e(t s - t K )\ + M K M V [e{t s ) + e(t s - t K )\ e(U - t.) + 2M 2 } , (8) 

where t s is the time-component of the space-time position of the strong chiral Lagrangian vertex in this 
diagram. The function e(t) is defined as: 

r +1, t > 0. 

e(t) = (9) 
I -1, t < 0. 

In the above two equations, only the integral of t s between and t„ can result in a "vanishing denomi- 
nator" when Mk — > Explicitly it is: 

Ml 



2(b)|0-»t, 



When Mr —* M^, the factor 
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(10) 



3 -2i(M K -A/^)t, 



-2i(M K - M w ) 

is just t % . Therefore, for finite t n (or finite in Euclidean space), C2 (or C2 in Euclidean space) vanishes 
at Mk = M^ (with both C*2( a ) = and C2(&) = separately) due to the explicit factor of Ml — Ml, and 
there is no power divergence. This conclusion remains true to all orders in xPT. 

To conclude, wc would like to discuss in some more detail why the factor linear in t n appears in Eq. 
QlUfl. even though C 2 (b) vanishes for Mk = M w because of the explicit factor (m s — m a ) in Eq. (J5J. 
Omitting this factor, our result contains a term linear in t n for Mk — M n . One would expect that if one 
takes large after taking the limit Mk — > M ff , it would be necessary to unitarize C 2 ( ). 2 Re-interpreting 
diagraml^b) as the lowest-order contribution in xPT (in the strong vertex) to K°K° — ► tt + it~ scattering 
(as we did above), the term linear in t„ can be understood as follows. For Mk = M„, there is full SU(3) 
symmetry, and \KK) and |7T7r) s-wave, 1 = states can be expressed in terms of the (1 = components 
of the) irreducible states |1), |8) and \27) through the relations ||] 



ID = ±\KK) 
|8) = ^\KK) + ±\ m )-^\w) 



|27) 



KK 




7T7T) 



(11) 



2 The on-shell divergence discussed earlier has nothing to do with this term linear in t^, but, as explained above, with 
the oscillatory behavior of e~ lJUff I'™ ~' s I . 
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where 



\KK) = ^=(\K°K°) + \K+K-)), 
Itttt) = -L (|ttV) + ^|tt+0) • (12) 
From these relations, it follows that 

(nn(t = U)\KK(t = 0)) = = t n )\l(t = 0)) - ^-(8(t = t„)\8(t = 0)) - = t,)\27(t = 0)). 

4 5 20 

. (13) 

To leading order in xPT this expression contains a term linear in the coefficient of which is the 
corresponding linear combination of finite- volume two-particle energy shifts, thus explaining how a term 
linear in t % appears in C 2 (b)- Note that the normalization of C 2 {b) has been chosen such that (in leading 
non-vanishing order) it is independent of the spatial volume. Higher orders in xPT will indeed unitarize 
our result. 
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